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I. INTRODUCTION 


Ball bearings are among the most important and 

frequently used components of the machines; however 
bearings may contain manufacturing errors or mounting 
defects. Such errors cause vibration, noise, and even failure 
of the whole system, which leads to expensive claims for 
damage. To avoid this and to ensure rapid and cheap 
production there is a need for a quick end-test of the 
machines to determine any bearing faults. The first step in 
bearing-fault detection during run-up would be a numerical 
model for the bearing-vibration response due to faults 
during the run-up. A well-defined vibration signal during 
the run-up of a faulty bearing could be used to find a 
suitable method for the fault diagnostic. A lot of research 
work has been done to model the vibration response of a 
bearing due to faults at a constant rotational speed. 
Tadina [ 1] presented an improved 2D bearing model for 
investigation of the vibrations of a ball-bearing during run- 
up. They presented numerical model assumes deformable 
outer race, which is modeled with finite elements, 
centrifugal load effects and radial clearance. 

Various surface defects due to local deformations 
are introduced into the developed model. The detailed 
geometry of the local defects is modeled as an impressed 
ellipsoid on the races and as a attended sphere for the 
rolling balls. The obtained equations of motion were solved 
numerically with a modified New mark time-integration 
method for the increasing rotational frequency of the shaft. 
The simulated vibration response of the bearing with 
different local faults was used to test the suitability of the 
continuous wavelet transformation for the bearing fault 
identification and _ classification. The experimental 
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determine the machine vibration response. The rolling bearing, with outer ring 
fixed, is a multi body mechanical system with rolling elements that transmit motion 
and load from the inner raceway to the outer raceway. The specific construction of 
a bearing has a decisive influence on its dynamic behavior. The paper defines a 
vibration model of a rigid rotor supported by rolling element bearings. By 
application of the defined model, the parametric analysis of the influence of 
internal radial clearance value and number of rolling elements influence on the 
rigid rotor vibrations in unloaded rolling element bearing was presented. 


Index Terms : analytical model, rolling bearing, ball passage frequency, 
internal radial clearance. 


validation of the simulation was done by wavelet 
transformation. 

Upadhyay [2] presents a mathematical model to 
investigate the nonlinear dynamic behavior of a high speed 
rotor-bearing system due to defects of rolling elements. 
They used two defects for the study which are off size 
rolling element and waviness of rolling element. In the 
formulation, the contacts between rolling elements and 
inner/outer races are considered as nonlinear springs and 
also used nonlinear damping, which is developed by 
correlating the contact damping force with the equivalent 
contact stiffness and contact deformation rate. The 
equations of motion are formulated using Lagrange’s 
equation, considering the vibration characteristics of the 
individual components such as inner race, outer race, 
rolling elements and rotor. They solved the equation of 
motion using numerical line integration technique. 

Rubio [3] developed rolling bearing analytical 
formulation, the contact between rolling element and 
raceways is considered as nonlinear springs and their 
stiffness are obtained by using Hertzian elastic contact 
deformation theory. Equation of motion is developed using 
Lagrange’s formulation and simulation is validated using 
Algor© code and finite element method in time domain 
with dynamic events. 

Wensing [ 4] considered the stiffness and damping 
of the elastohydro dynamic lubricated (EHL) contact 
between ball and guiding rings. They formulated the 
equation using finite element modeling and time integration 
and component mode synthesis (CMS) is used as method of 
solution of equation. Experimental validation is done by 
vibration test spindles. 

Tillema [5] presented the 3D nonlinear time 
dependent computational model of rolling element bearing. 
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They studied the vibration generation and _ vibration 
transmission characteristics of the rolling element. The 
formulation was done by Maxwell model of time domain 
equation. The equation of motion was solved by Newmark 
time integration method and Newton-Rapson method. The 
mathematical simulation was validated by experiments on a 
rotor dynamic test rig. 

Jang [6] used the two type of defect of the bearing 
to study the dynamic behavior of the bearing. Those were 
wavieness of ball and waviness of the races. Problem 
formulation was done by Runge-Kutta-Fehlberg algorithm 
and validation was done with those of prior researchers. 

Purohit [7] developed the nonlinear dynamic 
model of ball bearing with effect of varying preload and 
number of balls. The mathematical model was solved by 
Newmark-f with Newton-Rapson method. The results are 
validated using FFT. 

Harsha [8] presented nonlinear dynamic analysis 
of a rail axle-bearing system due to the effects of the 
number of rolling elements has been studied. The result 
showed The results show the appearance of instability and 
chaos in the dynamic response as the speed of the axle- 
bearing system is changed. Period doubling and mechanism 
of intermittency have been observed which lead to chaos. 
The appearance of regions of periodic, sub-harmonic and 
chaotic behavior is seen to be strongly dependent on the 
radial clearance. Poincaré maps, phase plots and frequency 
spectra are used to elucidate and to illustrate the diversity of 
the system behavior. 


Il. DEFECTS CONSIDER FOR STUDY 


The dynamic performance of bearing is highly 
influential on performance of any machine. More 
specifically, the presence of bearing defects often results in 
reduced efficiency, or even severe damage of the machine 
under consideration. System rigidity in machine tool is 
extremely important because the magnitude of deflection 
under load determines machining accuracy. The bearing 
stiffness is main factor that influence system rigidity. 
Therefore, bearing preload is needed to enhance system 
rigidity and to increase running accuracy. Bearing preload 
can be regarded as negative internal clearance. A proper 
amount of negative bearing clearance is desirable in order 
to stiffen the support of the spindle. However, inappropriate 
negative bearing clearance can cause excessive rolling 
contact stresses and eventually lead to bearing seizure. 
Therefore, Proper internal clearance have to be selected in 
order to prevent bearing seizure and to improve bearing 
stiffness. 


IN. THE EQUATION OF MOTION 
The complex motion of balanced rigid rotor sleeve 


in rolling element bearings and found that the trajectory of 
rotor cross-section center in bearings with the internal 
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radial clearance and idealized geometry could be 
represented by the equation of absolute sinusoid: 


y=A [sin “| (1) 
Where A is the peak-to-peak (pp) rotor vibration amplitude, 
i the arc distance between the bearing rolling elements and 
s the curviliniear coordinate of bearing rolling elements 
brunt. 

The path of rolling elements brunt is equal to pitch diameter 
of the cage, ie, centers of rolling elements are moving on 
the cage pitch circle. Based on this , the arc distance 
between the rolling bearing elements (A) can be obtained 
according to the following equation: 
tee 


A= (2) 


Zz 
Also, the curvilinear coordinate of bearing rolling elements 
brunt (s) can be expressed as 


1 
S=T.9 = 54.P (3) 


If the expression (3) and the expression (2), are applied to 
(1), after rearrangement: 


y = A|sin? o| (4) 
Angular coordinate of rolling elements brunt (@) can be 
expressed in function of time over cage frequency @, 
according to the following equation: 


p= at (5) 


But w, = “(4 = cosa), so 
g = at =2(1-Pcosa)t (6) 


Where w is shaft frequency, d, the diameter of 
rolling elements, a the bearing contact angle and t the 
time. 

hy 








Fig. 1. Schematic view of rolling element bearings kinematics. 
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Now the equation for curvilinear coordinated s can be 
given as 


s=+u,d,t =*(1-“cosa) (7) 


Furthermore, if the terms for s and A are added to eq. 
(1)one can obtained: 


y=A [sin (Ew.t)| =A [sin z2(4 — ‘cos a) || (8) 
Or 

y =4|sin()| (9) 
Where 


tax ={0(1-“cosa) (10) 

In above equationf,,,, as previously noted, represent the 
ball passage frequency (BPF). 

According to this, the rotor cross-section centre in relation 
to the centre opening in the bearing housing vibrates 
constantly during bearing operation. These vibrations can 
be represented by the equation of absolute sinusoid (9).The 
frequency of these vibrations is equal to ball passage 
frequency, and its amplitude is equal to rotor vibrations pp- 
amplitude. The pp-amplitude (A) is equal to the difference 
of maximum divergence from a reference rotor position on 
one and the other side 


A = Vmax — Ymin (11) 


IV. THEORETICAL ANALYSIS OF RIGID ROTOR 
VIBRATIONS AMPLITUDE CHANGE IN UNLOADED 
ROLLING ELEMENT BEARINGS 

In an unloaded rolling element bearing ,with 
internal radial clearance, the inner race is supported by two 
rolling elements as shown in Fig. 2. This support case often 
occurs in bearings subjected to low values of external radial 
load, in which the contact deformations are not sufficient 
enough to enable the participation of other rolling elements 
in the transfer of load. In such cases, the bearing structure 
can be observed as a system of rigid bodies ,in which the 
transfer of load from the inner to the outer race goes over 
one or two rolling elements. In a system of rigid bodies , all 
links are geometric and all movements can be analyzed 
based on geometrical laws. The rotor centre will vibrate 
between two extreme positions, which are determined by 
position of the set of rolling elements, as is shown in Fig. 2. 
As shown in Fig. 2a, sleeve centre will be moved in relation 
to the outerrace centre for half the value of internal radial 
clearance (e). This is the minimum deviation of the rotor in 
relation to the opening centre in the bearing case, i.e.: 
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Inn = Fate He a (12) 


The maximum deviation of the rotor in relation to the 
centre of the inner race occurs in a case of symmetrically 
distributed rolling elements in relation to the bearing 
vertical axis (Fig. 2b). The maximum value of rotor 
divergence Y,,,,, can be obtained based on the scheme 
shown in Fig. 3[28]. 

From Fig. 3 one can see that the following equality applies: 





aan 2 
OB, - OA? = [(*2*) — yz, sin? 4 (13) 


2 
2V be i d, bi 2 i 
> can be ignored, because Yngx 18 the 


BoA = 
Article y2,, sin 


small dimension of higher order in relation to the d; and dy, 
so 








(14) 





aA V 
OA=y._,. COS—; 
Pee 2 
O'B, _4, +d, é 
2 
0B, = =O +e; 








; _ 7 
O'A=y_. sin; 
sy enae 2 


Fig.3. Scheme of the maximum divergence of the rotor leaned on a two 
rolling bodies. 





If this is inserted in the equation OB, = OA + B,A: 





dj—d,+e 
2 


dj+dp 


: (15) 


y 
= Vmax COS > + 


After systemization the maximum deviation of the rotor 


Vmax : 
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e 1 

Vmax = 2 bos ("/>) (16) 

The difference of the extreme values of rotor divergence 
gives the maximum gap ,i.e .pp-amplitude (4,) of rigid 
rotor centre vibrations in the case of support by two rolling 
elements, according to 





1 
A, = Vmax ~ Vmin = <( yo 1) (17) 


cos > 
2 


By analyzing expression (17) one can see that the pp- 
amplitude of rigid rotor vibrations in rolling element 
bearings, in unloaded bearing, depends on the size of the 
internal radial clearance and angle distance between the 
rolling elements. Since the angle between rolling elements, 
is directly dependent on their number, one can say that the 
rotor vibration pp-amplitude, is directly determined by the 
value of the internal radial clearance and the total number 
of rolling elements in the bearing. Other bearing 
characteristics do not have influence on the size of pp- 
amplitude. If the external load is relatively low, bearings 
with the same number of rolling elements and the same 
value of internal radial clearance will have the same value 
of vibration, regardless of the size and dimensions of the 
bearings structure elements. 


Fig. 4. Dependence of k,, on d,/ d, ratio and number of rolling elements 
(z) for single row radial contact ball bearing. 


V. PARAMETRIC ANALYSIS OF THE INFLUENCE OF 
BEARING CONSTRUCTION ON BALANCED RIGID 
ROTOR VIBRATION CHARACTERISTICS 

A. Vibration Frequency 


Based on previous analysis it can be concluded 
that vibration frequency of rigid rotor in idealistic rolling 
element bearing, with internal radial clearance, is equal to 
ball passage frequency. The frequency change is given by 
Eq. (10) and there from it is clear that ball passage 
frequency depends directly on rotor frequency and bearing 
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construction (number of rolling elements z, ratio of rolling 
element diameter db and cage diameter dc and contact 
angle a). Furthermore, from Eq. (10) it is obvious that BPF 
does not depend on the value of internal radial clearance. If, 
in Eq.(10),use that 


ke =i (1-2) cosa (18) 


it can be seen that this frequency is linearly dependable on 
the shaft frequency, with linearity coefficient k,, 


fo = kaw (19) 


Coefficient k,, represent the linearity coefficient whose 
value is a function of bearing construction. This coefficient 
is a ratio of BPF and shaft frequency and shows how much 
BPF is larger than shaft frequency. k,, is a function of total 
number of rolling elements (z), ratio of rolling elements 
diameter and cage diameter (db/dc), as well as contact 
angle (a). By multiplication of k,, with shaft frequency, it 
is easy to obtain the BPF. Fig. 4 shows the dependence of 
coefficient k,, on the bearing construction for radial contact 
ball bearing. The diagram enables the obtaining of k,, value 
in respect to ratio db/dc, for different values of total number 
of rolling elements (z). By analyzing Fig. 4 it is clear that 
k,, linearly decreases with increase of ratio db/dc, while its 
value increases with the increase of total number of rolling 
elements. Based on bearing construction (number of rolling 
elements and db/dc ratio), the BPF can be as much as ten 
times larger than the shaft frequency. Fig. 5 gives the 
dependence of coefficientk,,, thus BPF, on the total number 
of rolling elements in radial contact rolling bearing. Based 
on Fig. 5 and Eq. (10) it can be concluded that BPF is linear 
function of total number of rolling elements, with linearity 
coefficient that depends on ratio db/dc and shaft frequency. 
Fig. 5 also shows that with the increase of number of 
rolling elements, k,, also increases and thus BPF. Steeper 
lines have lower values of (db/dc), i.e .the increase of 
number of rolling elements with smaller values of db/dc has 
greater influence on increase of BPF. 





z 
Fig. 5. Dependence of k,, on the number of rolling elements for different 
values of ratio between rolling element diameter and cage pitch diameter 
(db/dc). 
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Fig. 6. Three-dimensional diagram of BPF dependency in relation to value 
of internal radial clearance and total number of rolling elements in the case 
of reclining on two rolling elements. 


B. Vibrations Amplitude 

The influence of bearing construction on BPV amplitude, in 
unloaded bearing, can be analyzed by Eq. (17) which 
defines the pp-amplitude of balanced rigid rotor centre 


vibrations in a case of support by 
25 


20 





0 20 40 60 80 100 120 140 160 180 200 
e, um 


Fig.7. Balanced rigid rotor vibration amplitude dependency in relation to 
value of internal radial clearance in the case of support on two rolling 
elements 


two rolling elements. According to (17), amplitude linearly 
depends on the value of internal radial clearance, with 
linearity coefficient that is a function of angle between 
rolling elements (y): 





ae 1) (20) 


1 

os tg(a) = (<> 

where a is a linearity coefficient. The total number of 
elements in the bearing determines the angle between 
rolling elements (y), so the linearity coefficient is direct 
function of the total number of rolling elements. Therefore, 
the analysis of the internal radial clearance impact on BPV 
amplitude is best to perform for different values of total 
number of rolling elements. Fig. 6 gives the three- 
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dimensional dependency of pp- amplitude in relation to 
value of internal radial clearance and total number of 
rolling elements. Fig. 7 gives the dependency of pp- 
amplitude on the internal radial clearance for different 
values of total number of rolling elements. With the 
increase of value of internal radial clearance ,the value of 
pp-amplitude also linearly increases. The increment is 
bigger with the decrease of total number of rolling 
elements. From Fig. 7 it is clear that with the increase of 
total number of rolling elements the dependency lines tend 
towards horizontal, i.e. the angle that they precede in 
relation to horizontal axis decreases. 


VI. CONCLUSIONS 


Based on the numerical analysis it can be concluded that: 

1. The centre of rotor cross-section continuously vibrates 
in regard to bearing housing opening. These vibrations 
can be represented by the equation of the absolute 
sinusoid. By this equation, the vibration frequency of the 
rigid rotor in the ideal rolling element bearing with 
internal radial clearance is equal to ball passage 
frequency (BPF). The amplitude of these vibrations is 
equal to rotor vibrations pp-amplitude and it is equal to 
the difference of maximum divergence from a reference 
rotor position in one and the other side. 

2. The BPF is in linear relation to the shaft frequency, with 
the coefficient of linearityk,,. Coefficient k,, shows the 
impact of the rolling element bearing construction 
parameters on the BPF. This coefficient depends on the 
following: number of rolling elements (z), ratio of rolling 
elements diameter and the cage pitch diameter (db/dc), as 
well as the contact angle of the bearing (a). The value of 
internal radial clearance does not influence k,, 

3. The BPF linearly increases with the increase of number 
of rolling elements. This increase is more pronounced 
with the lower values of ratio between diameter of rolling 
elements and cage pitch diameter (db/dc). 

4.The vibration amplitude on BPF, in unloaded rolling 
element bearing ,is linearly dependable on the value of 
internal radial clearance, with linearity coefficient, which 
is influenced by angular distance between rolling 
elements, i.e .the total number of rolling elements in the 
bearing. With the increase of the internal radial clearance, 
the value of amplitude increases linearly. The increase 
gets much bigger as the total number of rolling elements 
decreases. 


NOMENCLATURE 


trajectory of rotor cross-section centre 

peak-to-peak (pp) rotor vibration amplitude 

arc distance between the bearing rolling elements 
curvilinear coordinate of bearing rolling elements brunt 
angular coordinate of rolling elementsbrunt 
total number of rolling elements 


Neg 2% >boX< 
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time 

shaft frequency 

cage frequency 

contact angle 

ball passage frequency (BPF) 
calculated values of BPF 


fope eXperimentally obtained values of BPF 


y angle between rolling elements 

e internal radial clearance 

d nominal diameter (shaft diameter) 

di bearing inner raceway diameter 

dc cage pitch diameter 

db diameter of rolling elements 

do bearing outer raceway diameter 

D bearing outside diameter 

ri inner raceway groove curvature radius 

ro outer raceway groove curvature radius 

kw coefficient a coefficient 

a coefficient 
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